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ABSTRACT 

We study the properties of two bars formed in fully cosmological hydrodynamical 
simulations of the formation of Milky Way- mass galaxies. In one case, the bar formed 
in a system with disc, bulge and halo components and is relatively strong and long, as 
could be expected for a system where the spheroid strongly influences the evolution. 
The second bar is less strong, shorter, and formed in a galaxy with no significant 
bulge component. We study the strength and length of the bars, the stellar density 
profiles along and across the bars and the velocity fields in the bar region. We compare 
them with the results of dynamical (idealised) simulations and with observations, and 
find, in general, a good agreement, although we detect some important differences as 
well. Our results show that more or less realistic bars can form naturally in a ACDM 
cosmology, and open up the possibility to study the bar formation process in a more 
consistent way than previously done, since the host galaxies grow, accrete matter and 
significantly evolve during the formation and evolution of the bar. 

Key words: galaxies: kinematics and dynamics - galaxies: structure - cosmology: 
theory - galaxies: photometry - methods: numerical 



1 INTRODUCTION 

Bars are present in roug h ly two thirds of a ll galactic 
discs (jEskridge et al.1 l200ol : iBarazza et all l2008h and can 
drive the secular evolution of their host galaxies. They 
have thus been the subject of a number of studies based 
on dynamical simulations, i.e. simulations aimed towards 
an understan ding of the main re le vant dynamical mech- 
anisms (e.g. |Combes et al. 1990l: iDebattista &; Sellwoodl 



|2000l: IXt hanassoula & Misi riotisl 120021 : lAthanassoulal 2003: 
iMartinez-Vaipuesta et al.l 12006). These use idealised initial 
conditions, made specifically for the question under study, 
and set aside all other effects, at least until an understand- 
ing of the basic mechanism is achieved. The initial conditions 
correspond to a fully developed disc plus halo system which 
is as near equilibrium as possible, so that the bar formation 
can be studied uninfluenced by other instabilities. Further- 
more, interactions with other galaxies, major or minor, as 
well as inflow from the en vironment has bee n seldom taken 
into account (see, however. I Curir et alj|200q ). while the halo 
is usually assumed to be initially spherical. Moreover, the 
gas is generally neglected or modelled without taking into 
account star formation, feedback or cooling. 

In this letter, we take a different approach, i.e. we study 
bar formation in the context of the ACDM cosmology. In 
this way, we include important external effects, such as ac- 
cretion and interactions and, more important, the bars do 



not wait for the galaxy to be fully formed to start their own 
formation. The disc and even the dark matter halo keep 
growing while the bar forms. Also, our galaxies do not have 
preset halo-to-disc mass ratios, velocity dispersions, or other 
properties, and all their properties directly result from the 
simulations. This of course implies that we will not be able 
to address the same questions as the dynamical simulations. 
For example, we can not examine how a given property of 
the host galaxy, e.g. the halo radial density profile, will influ- 
ence the formation and evolution of the bars. On the other 
side, we will include in our simulations more physics that 
in any other single bar formation simulation. Our main goal 
is to investigate if bars can naturally form in ACDM, and 
compare their properties with those of dynamic simulations 
and with observations. 

The layout of the paper is as follows. In Section [2] we 
describe the simulations, in Section [3] we present and discuss 
our results, and we conclude in Section [4] 



2 SIMULATIONS AND ANALYSIS 

We study the present-day properties of bars formed in two 
hydrodynamical simulations of galaxy formation in a ACDM 
universe. The simulations correspond to galaxies that, at 
redshift z — 0, are similar in mass to the Milky Way and 
are mildly isolated. They are a su b-sample of the eig h t sim- 
ulations extensively described in IScannapieco et aL {2009, 
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Figure 1. Projected surface stellar density for Aq-C (left-hand 
panels) and Aq-G (right-hand panels). From top to bottom, the 
figures show the face-on (XY), side-on (XZ) and end-on (YZ) 
views. The black lines correspond to isodensity contours equally- 
spaced in logX, and the white circles indicate the bar length ob- 
tained with three different methods (see Section [3~4|) . 



1201(1 l201ll . hereafter S09, S10 and Sll, respectively). We 
analyse here the galaxies named Aq-C (Aq-C-5 in S09, 
M200 = 1.6 x 1O 12 M ) and Aq-G (Aq-G-5 in S09, M 200 = 
6.8 x 10 n M©), where bars are clearly present at z — 0. Two 
other galaxies, Aq-A and Aq-E, have bars at z = (S10); 
however, we defer their analysis to a separate work since the 
evolution of these systems is far more complicated. Aq-A has 
two stellar misaligned discs at z = 0, while Aq-E has a ro- 
tating bulge and a recent interaction with a satellite galaxy 
which strongly disturbs the system (Sll). 

We assume a ACDM cosmology with Qa — 0.75, Q m — 
0.25, Q b = 0.04, a 8 = 0.9 and H = 73 km s" 1 Mpc" 1 . 
The mass resolution is ~ 10 6 M© for dark matter and ~ 
2 x 1O 5 M0 for baryons (S09). The gravitational softening is 
1.4 kpc, it is fixed in comoving coordinates and is the same 
for gas, stars and dark matter particles. 

The sim ulations were run w ith the Tree-PM SPH code 
GADG ET-3 (ISpringel et al. | l2008l) with the additional mod- 
ules of lScannapieco et al.1 120051 12006). It includes star for- 
mation, chemical enrichment and (type II and la) supernova 
feedback, metal-dependent cooling, an explicit multiphase 
model for the gas component, and the effects of a UV back- 
ground. We refer the reader to S09 and Sll for full details on 
the initial conditions, simulation code and resolution effects. 

Our analysis meth ods largely follow those 
lAthanassoula k Misiriotid (|2QQ 



easily compare our results to dynamic simulations and to 
observations. There are, however, two differences. First, 
unlike in dynamic simulations where it is possible to analyse 
discs separately from bulges (all through this paper, by 
bulge we mean classical bulge, iKormendv &; Kennicuttl 
2004), here they constitute together the stellar compo- 
nent and any attempt to distinguish between them is 
approximate. The second difference is that in dynamic 
simulations the softening is about one tenth of that used 
here. Consequently, our results will be more smoothed out 
than those of dynamical simulations. 

We will also follow the terminology introduced in AM02, 
where MH denotes simulations with strong bars, in which 
the near-resonant material in the bar region emits a con- 
siderable amount of angular momentum, which is absorbed 
mainly by the near-resonant material in the halo. Simu- 
lations in which considerably less angular momentum has 
been redistributed, a nd which therefore have less strong bars 
(|Athanasso ula 2003, hereafter A03), are called MD. 

3 RESULTS 
3.1 Morphology 

Fig. [1] shows maps of projected surface stellar mass density 
for Aq-C and Aq-G. The projection is such that XY is the 
disc's equatorial plane, with angular momentum in the pos- 
itive Z direction and with the bar major, intermediate and 
minor axes lying along the X, Y and Z axes, respectively. 

Aq-C has a significant well-defined disc, a bulge-like 
component and a bar. The bar looks very strong, very thin 
and rectangular-like. Seen face-on, the bulge has a rather pe- 
culiar shape. In contrast, Aq-G does not have a significant 
bulge (S10), but has an important bar and a clearly iden- 
tifiable disc. Its bar is very symmetric and, viewed face-on, 
much fatter than that of Aq-C. 

We find some differences between our simulated galax- 
ies and both real galaxies and galaxies in dynamic simu- 
lations, in particular in the face-on views. We find that 
Aq-C has a bizarre shaped bulge, unlike what is usually 
observed. Note however that, in a ACDM cosmology, the 
shape of galaxies can change very rapidly particularly during 
and/or after interactions with satellites and mergers. In the 
case of Aq-G, the disagreement is quantitative, rather than 
qualitative. Rectangular-shaped bars ha ve been many times 
observed, but always in strong bars ( Athanassoula et al.l 
Il990l : lAthanassoula fe Misirioti3 120021 : badottil l201lh . On 
the other hand, the bar in Aq-G is fat, more like an oval 
(i.e. weak), but still has a clear rectangular outline. Aq-G 
then looks like a hybrid between strong and weak bars. 



3.2 Bar strength 

In order to quantify the bar strength we Fourier analysed 
the projected face-on mass density (£), i.e.: 



E(i*,0) 



A (R) 



+ (i^)cos(m<9) + ^ (i^)sin(m(9)] (1) 



with being the azimuthal angle and R the cylindrical ra- 
dius. In practice, we calculate A m (R) and B m (R) as: 



hereafter AM02), 



m 

to 



A m (R) — ^2 m% cos ( m 6i)i 171 > 



(2) 
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Figure 2. Le/fc: Relative Fourier amplitudes of the even compo- 
nents up to m = 8. Right: Projected mass density profiles along 
the bar major axis (X, dashed lines), the bar minor axis (Y, dot- 
ted lines), and azimuthally averaged (solid lines), when the galaxy 
is seen face-on. 



B, 



.(«) = E 



i(mOi), m > 



(3) 



where m l is the mass of the stellar particle i and the sum- 
mation is over all particles in an annulus around radius 
R. We define the amplitudes of the Fourier components as 
Ira — \J + B% } (m > 0) and io = Ao/2, and use their 
ratio as a measure of the bar strength (note that several 
dynamical simulation studies and observations use the ra- 
tios Im/Ao instead of Im/Io.). We show Im/Io for Aq-C and 
Aq-G (up to m = 8) in Fig. [2] (left-hand panels). 

Aq-C, the strongest of our two bars, has a large 
m = 2 component. Its maximum, (i2//o) max — 0.92, oc- 
curs at i?max = 4.5 kpc and is similar to those found 
in dynamical simulation s with gas but no star formation 
(| Villa- Vargas et aT1l2010l and references therein). Compared 
to dynamical simulations with star formation, feedback and 
cooling, Aq-C corresponds roughly to simulations with most 
baryons initially in a gaseous component and reaching about 
6 to 10% gas in the disc at times comparable to z = 
(Athanassoula et al., in prep.). The shape of the radial pro- 
file is more reminiscent of that found in gas-rich dynamical 
simulations. Some of the differences could, nevertheless, be 
due to the fact that we are including in the Fourier analysis 
all stars, i.e. from the disc, the bar and the bulge. This will 
necessarily lower the m components, though, in this case 
at least, not all m in the same way since the bulge has a 
clear rect angularity. Also, m — 4 aside, this lowering must 
be stronger in the inner parts, where the bulge contribution 
is larger. 

The higher order even moments (m = 4, 6, 8) have max- 
ima of h/Io = 0.64, J 6 /Jo = 0.44 and I 8 /I = 0.28, in very 
good agreement with MH-type simulations. These maxima 
occur at radii roughly equal to that of the m = 2, while 
in MH-type dynamical simulations they occur considerably 
further out. This also could be due to the strong rectangu- 
larity of the bulge (Fig. [1]). 

The bar of Aq-G is less strong, with (/ 2 //o) max = 0.70 
that occurs at i? ma x = 3.4 kpc. The m = 4, 6, 8 moments 
are very small (< 0.18) at all radii, in good agreement with 
the more elliptical- like bar outline. In fact, the m > 4 com- 
ponents are within the noise. Thus they are in many ways 
similar to those from MD simulations, where less angular 



momentum has been redistributed within the galaxy than 
in MH- types (AM02, A03). The fact that i? max is shorter 
in Aq-G than in Aq-C indicates that the bar is shorter, as 
already inferred from Fig. [TJ Note, however, that the two 
galaxies have neither the same mass, nor the same extent, 
as seen in Fig. [1] (see also Section f3.4p . 

Another important difference between the two bars 
is the shape of the I2/I0 amplitude for R < i? ma x and 
R > i?max- In Aq-C, the declines on the two sides of i? ma x 
are similar, whereas in Aq-G the decline is clearly steeper 
inwards. In dynamical simulations, the norm is an asym- 
metrical decline on the two sides of the maxim um (but see 
simulations MH1 in I Athanassoula et al.l [20051 , which has a 
strong central mass concentration). 

Finally, the odd Fourier components (m = 1, 3, 5, 7) are, 
both for Aq-C and for Aq-G, very small (< 0.1) at all radii 
indicating that both bars are quite symmetric. 

Our results are in relatively good agreement with ob- 
servations, that find a great variety of I2/I0 profi les and of 
maximum values, typically between 0.4 and 0.8 ( Ohta et "all 
Il990l : lAguerri et alll200lL Bool iButa et al J 120061 ). In AchC. 
/ 2 //™ ax = 0.92 is v ery high, but val ues as high as 0.9 have 
also been observed (|Buta et alJl2Q06h . 



3.3 Density profiles along major and minor axes 

Fig. [2] (right panels) shows the stellar density profiles viewed 
face-on, along the major (X) and minor (Y) axes of the bar, 
as well as azimuthally averaged (note that we include here 
all stars in the simulations, i.e. bulge, disc and bar). 

The profiles of the two simulations differ significantly 
between them, while having a number of similarities with 
those of dynamic simulations. The profile of Aq-G along 
the bar major axis shows a clear near-exponential drop; 
a behaviour similar to that of MD models, and to the 
' exponential' bars often obs erved in late-type galaxies 
(Elm egreen &; Elmegreenll 19851 ). In contrast, Aq-C has a less 
typical profile. It shows a flattish section, but of short ex- 
tent, arguing that it is linked to the bulge rather than the 
bar, contrary to the flattish parts of MH bars (Fig. 5 in 
AM02) and of 'flat' bars in observed early-type galaxies 
([Elmegreen &; Elmegreen! [l985h . Furthermore, and contrary 
to MH-type bars in dynamic simulations, there is no abrupt 
fall of the density profile at the position of the bar end, 
presumably due to the fact that the softening in the cos- 
mological simulations is much larger and so washes out the 
relatively steep drop at the end of the bar. 

The density profiles along the minor axis of the two 
galaxies also show significant differences. In Aq-C, the profile 
shows a clear change in slope at the position where the bar 
ends. This behaviour is also found in MH-type bars and 
means that the bar is sufficiently strong to clear out the 
region around the Lagrangian points, as one sees in many 
observed strong barred galaxies. Aq-G also shows a change 
in the slope of the density profile along the bar minor axis; 
but much less pronounced than Aq-C. Fig|2]also shows that 
the difference between the cuts along the major and minor 
axes in the bar region is considerable larger in Aq-C than 
in Aq-G, confirming that the bar is much stronger in the 
former case. 
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Figure 3. Phase of the stellar distribution as a function of radius. 
The dashed lines indicate the phase of the bar (</>b), and the 
dotted lines are at </>b± arcsin(0.3). 



3.4 Bar length 

AM02 presented a number of ways to measure bar length 
and discussed their advantages and disadvantages. When- 
ever possible or useful, these me thods were extended to ob- 
served bars (|Gadotti et al.| [2QQ7). As not all of these meth- 
ods are well adapted to cosmological simulations, we will 
here use three of them, of which one after a small modifica- 
tion. 

The first method (Hi in AM02, bar length L p hase) uses 
the radial profile of the phase of the m — 2 component 
which, in the bar region, should be approximately constant. 
Outside it, the phase either shows a coherent increase or 
decrease due to a spiral, or large variations in the absence 
of clear structures. Following AM02, we calculate the phase 
of the m — 2 Fourier component of: (i) the whole bar-disc 
system and (ii) as a function of radius, and define the bar 
length as the maximum radius where these two quantities 
differ by less than arcsin(0.3). In this way, we obtain bar 
lengths of 8.8 kpc for Aq-C and 8.0 kpc for Aq-G (see Fig. [3] 
and dashed circles in Fig. [TJ. 

For the second method (iv in AM02, bar length L m =2) 
we use the relative m — 2 Fourier component. If the disc and 
bar were rigid, the length of the bar would be the radius at 
which the m — 2 component goes to zero. This is not the 
case in simulated bars, but we can still get an estimate of 
the bar length from the radius where the m — 2 amplitude 
drops to a given (arbitrary) fraction of the maximum, that 
we take to be 25%. In this way, we obtain lengths of 8.4 kpc 
for Aq-C and of 7.7 kpc for Aq-G (solid circles in Fig. [T]). 

The third method (v in AM02, bar length L pro f) uses 
the density profiles shown in Fig. [21 Namely, we take the 
difference between the profiles along the major and minor 
axes, which is zero at the centre, increases until a maximum 
and then drops again. If the disc and bar were rigid, the end 
of the bar would be where the two projected density profiles 
become equal again. In simulated bars, the bar length is de- 
fined as the radius where the difference between the profiles 
drops to 5 per cent of the maximum. This method yields 
bar lengths of 8.5 kpc for Aq-C and of 6.8 kpc for Aq-G 
(dot-dashed circles in Fig. [1]). 

The bar lengths obtained by our three methods agree 
very well between them for Aq-C and reasonably well for 
Aq-G. For Aq-C we get a mean value of ~ 8.5 kpc and for 
Aq-G of ~ 7.5 kpc. 

We can compare our results t o the ob s ervati ons of 
~ 300 barred galaxies presented in iGadottil (|201lL here- 
after Gil). By ap plying bulge-di sc-bar decompositions with 
the BUDDA code (Gadotti 2008), Gil obtained bar lengths 




Figure 4. Velocity fields of the stellar component for Aq-C and 
Aq-G, for three orientations of the bar: along the ordinate (upper 
panels) , at 45 degrees to it (middle panels) and along the abscissa 
(lower panels). The line of sight is always along the ordinate. The 
colors represent the mean velocities along the vertical axis and 
the white lines are the corresponding isodensity contours. 



(£bar) and ratios between the bar length and the disc scale- 
length (Lhar/h). We have applied the same method to our 
simulated galaxies (S10), obtaining the values of the disc 
scale lengths. We find Lb&r/h values of 0.80— 1.17 for Aq-C, 
and of 0.89— 1.15 for Aq-G (depending on the bar length es- 
timation) . These are at the low end of the distribution given 
in Gil, but still compatible with it. On the other hand, if we 
make comparisons directly with the values of bar length in 
kpc, we find that they are at the high end of the distribution 
in Gil but, again, still compatible with it. 

3.5 Kinematic properties 

The most telling characteristics of a bar in a radial velocity 
field are the presence of kinematical axes which are not per- 
pendicular to each other, as well as a zigzag (Z) behaviour of 
the isovelocity curves corresponding to the systemic veloc- 
ity of the galaxy. This Z aspect is best seen in velocity fields 
where the bar makes an angle of 45 degrees with the major 
axis of the galaxy (e.g. Athanassoula 1984). To investigate 
the kinematic structure of our simulated galaxies, we show 
in Fig. [4] their stellar velocity fields, obtained as in AM02, 
i.e. in a way that is the most convenient for comparing with 
observations of galaxies at intermediate inclinations, and 
also allows direct comparison with the AM02 results. For 
Fig. |4] we rotate the galaxy according to the desired view- 
ing angles, project all particles on the equatorial plane, and 
observe their line-of-sight (i.e. along the corresponding plot- 
ted ordinate axis) velocity component. For the three views 
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shown in the figure, this is equivalent to observing along the 
bar major axis, at 45 degrees to it and along the bar minor 
axis, respectively. Note that Fig. |4] shows the stellar velocity 
fields, which are not directly comparable to gaseous ones. 

The zigzag (Z) behaviour mentioned above is clearly 
seen in our two simulated galaxies. The effect in Aq-C is 
quite strong, and is comparable to that in model MH of 
AM02. The effect is weaker in Aq-G, but still stronger than 
the corresponding MD model of AM02. When we observe 
along the bar major axis, the isovelocities show the standard 
spider diagram form. The central concentration is less than 
in the dynamical simulations of AM02, but this is due to our 
larger softening. Finally, when we view along the bar minor 
axis, the velocity fields have an inner region which roughly 
delineates the bar and where the velocity is lower. This was 
also found in AM02 and is due to the fact that near the 
ends of the bar the particles are at their apocenters, i.e. the 
velocities are lower. The mean velocities in those regions can 
be further lowered if the c orresponding period ic orbits have 
loops at their apocenters (| Athanassoulal Il992h . 



4 CONCLUSIONS 

We presented the first study of bars formed in fully cos- 
mological, hydrodynamical simulations of Milky Way-sized 
haloes. The simulations include star formation, metal- 
dependent cooling, feedback from supernova and a UV back- 
ground field. In particular, we investigated the morphology, 
strength and length of the bars, their projected density pro- 
files and kinematical properties. 

The strongest of our two bars formed in a bulge- disc- 
halo system (Aq-C), while the weakest in a galaxy with a 
disc but no significant bulge (Aq-G). Compared to dynam- 
ical simulations, the strongest bar is similar to those found 
in systems where a considerable amount of angular momen- 
tum is exchanged between the bar and the halo. In contrast, 
the weakest bar is more reminiscent of dynamical simula- 
tions where considerably less angular momentum has been 
redistributed. The bar strength difference between our two 
simulated bars could be due to the effect of the bulge on the 
angular momentum exchange (which is negligible in Aq-G), 
as found in pure N-body dynamical simulations of isolated 
discs (AM02, A03), but the effect of interactions and/or of 
the gas component could also be decisive. 

The Aq-C (strongest) bar is very thin and rectangular- 
like, while that of Aq-G (weakest bar) is very fat, and 
also rectangular-like. Rectangular shapes are observed in 
strongly barred galaxies, so Aq-G is, in this respect, like an 
hybrid between strong and weak bars. Another difference is 
that the bulge of Aq-C has a rather peculiar shape which 
has so far not been observed. The peculiar shape is presum- 
ably a result of the cosmological formation of our galaxies, 
where mergers and/or interactions are common at all times. 

The lengths of our bars are in relatively good agreement 
with observations of Gil, although in the long tail of the 
distribution. The strongest bar, formed in our most massive 
galaxy, is also the longest. On the other hand, the ratios 
between bar length and disc scale-lenght are at the low end 
of the observed distribution, but still compatible with it. 

The density profiles are similar to those found in dy- 
namical simulations; however, we detect some differences 



particularly in Aq-C and in the very central regions, prob- 
ably because we include stars from all components in our 
analysis, unlike in dynamical simulations. Finally, the kine- 
matic properties of our bars are similar to those observed 
and to those found in dynamic simulations. 

The fact that the bars were obtained in a cosmological 
setting, and that their properties agree relatively well with 
known properties of bars, is not minor. Our galaxies grow 
significantly from z — 3 to z — 1, where the bars are first 
detected in the simulations. Moreover, external (accretion, 
interactions and mergers) and internal (cooling, star forma- 
tion and feedback) effects strongly affect the galaxies during 
their evolution. For these reasons, our results are quite sat- 
isfactory and very encouraging, and should incite more work 
on the study of bar formation in a cosmological context. 



REFERENCES 

Aguerri J. A. L., Beckman J. E., Prieto M., 2001, ApSS, 
276, 691 

Aguerri J. A. L., Debattista V. P., Corsini E. M., 2003, 

MNRAS, 338, 465 
Athanassoula E. 1992, MNRAS, 259, 328 
Athanassoula E. 2003, MNRAS, 341, 1179 (A03) 
Athanassoula E., Lambert J. O, Dehnen W., 2005, MN- 
RAS, 363, 496 

Athanassoula E., Misiriotis A., 2002, MNRAS, 330, 35 
(AM02) 

Athanassoula E., Morin S., Wozniak H., Puy D., Pierce 
M. J., Lombard J., Bosma A., 1990, MNRAS, 245, 130 

Barazza F. D., Jogee S., Marinova I., 2008, ApJ, 675, 1194 

Buta R., Laurikainen E., Salo H., Block D. L., Knapen 
J. H, 2006, AJ, 132, 1859 

Combes F., Debbasch F., Friedli D., Pfenniger D., 1990, 
A&A, 233, 82 

Curir A., Mazzei P., Murante G., 2006, A&A, 447, 453 
Debattista V. P., Sellwood J. A., 2000, ApJ, 543, 704 
Elmegreen B. G., Elmegreen D. M., 1985, ApJ, 288, 438 
Eskridge P. B., Frogel J. A., Pogge R. W. et al., 2000, AJ, 
119, 536 

Gadotti D. A. 2008, MNRAS, 384, 420 

Gadotti D. A. 2011, MNRAS, 415, 3308 (Gil) 

Gadotti D. A., Athanassoula E., Carrasco L., Bosma A., 

de Souza R. E., Recillas E., 2007, MNRAS, 381, 943 
Kormendy J., Kennicutt Jr. R. O, 2004, ARA&A, 42, 603 
Martinez- Valpuesta I., Shlosman I., Heller O, 2006, ApJ, 

637, 214 

Ohta K., Hamabe M., Wakamatsu K.-L, 1990, ApJ, 357, 
71 

Scannapieco O, Gadotti D. A., Jonsson P., White S. D. M., 

2010, MNRAS, 407, L41 (S10) 

Scannapieco O, Tissera P. B., White S. D. M., Springel V., 

2005, MNRAS, 364, 552 

Scannapieco O, Tissera P. B., White S. D. M., Springel V., 

2006, MNRAS, 371, 1125 

Scannapieco O, White S. D. M., Springel V., Tissera P. B., 

2009, MNRAS, 396, 696 (S09) 
Scannapieco O, White S. D. M., Springel V., Tissera P. B., 

2011, MNRAS, 417, 154 (Sll) 

Springel V., Wang J., Vogelsberger M. et al., 2008, MN- 
RAS, 391, 1685 



© 0000 RAS, MNRAS 000, 000-000 



6 Scannapieco & Athanassoula 

Villa- Vargas J., Shlosman L, Heller C, 2010, ApJ, 719, 
1470 



© 0000 RAS, MNRAS 000, 000-000 



